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[57] ABSTRACT 

A parallel digital computer with which real-time simu- 
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lations can be performed. A multiplicity of elemental 
computers are joined in a parallel architecture through 
a master computer and shared memory system. Real- 
time operation of the elemental computers and master 
computer are time interleaved without overlap. Input 
data entering each elemental computer is processed 
according to the designated function during one time 
segment, succeeded by a master computer operating 
time segment during which the resultant elemental 
computer data is selectively transferred to the appropri- 
ate elemental computers for the next processing se- 
quence. A rea]*time clock synchronizes the alternating 
time segments, allocating sufficient time for the process- 
ing performed by the elemental computer and the subse- 
quent data transfer steps. Integration, unlike addition, 
subtraction and the numerous other arithmetic and 
logic functions, is performed only once per frame time. 
To maintain simulation accuracy, the frame time is 
defmed to be short in comparison to the period of the 
highest frequency present in the real-time simulation. 
On the other hand, the frame time must be sufficiently 
long to satisfy the time demands of the intra-frame time 
cycling between elemental and master computers which 
is bounded by the elemental computer processing rate 
and the number of non-integrator elemental computer 
functions between successive elemental computers op- 
erating as integrators. 



12 Claims, 6 Drawing Figures 
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FIG. 2 
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of layout organization or conventional compiling of the 
PARALLEL DIGITAL COMPUTER input program. Both the ECs and the MC communicate 

ARCHITECTURE with the MC memory to store, retrieve, and selectively 

transfer data moving within the parallel digital com- 
RIGHTS OF THE GOVERNMENT 5 puter (PDC). 

The invention described herein may be manufactured The real-time clock alternately ^^^^^^ J^^^ 
and used by or for the Government of the United States group then the MC indivulually. TTiereby <»a U « se- 
for all governmental purposes without the payment of quentially processed, for instance, by havmg the MC 
any royalty. resMlts from the immediately preceeding EC 

10 run sequence and selectively transfer that data among 
CROSS REFERENCE TO RELATED the memory and EC elements in preparation for the 

APPLICATION next EC sequence. The cycling between the MC and 

This is a continuation of my appUcation Ser. No. ECs is performed at a rate significantly greater than the 
06,055.424 fded July 6, 1979. now abandoned. highest frequency present in the simultation being per- 

15 formed, so that the accuracy of the digital integration is 
BACKGROUND OF THE INVENTION degraded. Whenever individual ECs ar actuated by 

This invention relates to a parallel digital computer the RTC, they undertake their designated functions at a 
architecture and more particularly, to a computer for rate set by an independent, cycle clock having a fre- 
performing programs digitally but in the format of an quency orders of magnitude greater than the system 
analog computer. ^ RTC itself. 

Analog and hybrid computers have been the tradi- nPsmiPTlON OF THE DRAWINGS 

tional tools for performing real-time simulations of DESCRIPTION OF THE UKA wiinuj> 

body or structural dynamics. Nevertheless, many attri- fig. 1 is a schematic block diagram of the parallel 
butes of the digital computer make it an appealing sub- digital computer (PDC) as a complete system, 
stitute for the analog or hybrid forms, when the capabil- 25 j schematically shows the functional block 

itiesofthe digital computer and the nature of simulation structure of an elemental computer EC and contem- 
admit. The fundamental deficiency of the conventional plated variations of the registers, 
digital computer architecture, in the performance of pjQ 3 ^ schematic diagram depicting the various 
time related simulations, resides in the cycle time. Since operator control functions avaUable with an embodi- 
time dependent simulations in such digital computers 30 PUC. 

require repeated looping or cycUng through the pro- ^ schematically depicts the timing sequence of 

gram, the highest simulation frequency attainable is the real-time clock (RTC). 

severely constrained by both die complexity of the ^ ^ schematic of an analog simulation upon 

simulation and the cycle time of the digital , superimposed its parallel digital computer 

Furthermore, it is nearly unpossible to overcome this 35 ILJ^^ ^^,,^!Z^^Alt^,ntur^ 

denciency by dividing a single real-time problem into ^^.*?KSS^p~gram diagram of .he PDC 

subprograms of controUed computation time, because , l^tuJ 7i^J^*i^.sfT7ir^ « 
the varSble interval integration algorithms used vary in f<^^ ^he sunulation of FIG. 5. 
processing time from cycle to cycle. Thus, it may gencr- DETAILED DESCRIPTION 

ally be said that the art of real-time simulation in con- 40 « , % 1 „♦ , n>r%n\ ^^^uif^t^^r^ 

ventional architecture digital computer systems will be P^^^.^ ^^^"^ "^T^Th^H.^^^l^^^^ 

severely limited by the cycle sp^ of available and disclosed herem 

foreseeable digital computers. O^^ other hand, com- conventional digital computers by defining an architec- 
plex sunulations at other than real-time can be divided which pertmts concurrent, panOlel operations to 

into smaUer elements to reduce the nmning time of 45 und^e real-tiine simulations. Broadly stated, the 
conventional digital computers One such approach is PEM: architecture includes a plurahty of elemental com- 
taught in U.S. Pat. No. 3,346.851 issued to J. E. Thorn- Peters (ECs). typified by a state-of-the-art microproc^- 
ton et al. interconnected through a master computer (MC) 

and its memory, and synchronized to a real-rime clock 
BRIEF SUMMARY ^ (RTC). This general structure appears in the schematic 

The present invention is directed to an architecture block diagram of FIG. 1. As shown, master computer 
for digital computers, by which architecture the com- central processing unit (MC CPU) 1. MC memory 2 and 
puters are structurally oganired and time synchronized multiple ECs 3 communicate by means of data bus 4 and 
to undertake real-time simulations of the type previ- address bus 6. The actual distribution of daU registers 
ously attainable only with analog or hybrid computers. 35 between MC memory 2 and the individual ECs can be 
In one form, the structural organizations consists of a lumped in the MC memory or distributed among the 
master computer (MQ. a MC memory, a multiplicity of ECs; in either case they are treated as being shared, 
elemental computers (ECs). and a real-time clock Real-time clock 7, present in the figure, synchronizes 
(RTC) synchronizing in alternating fashion the EC the operations of MC CPU 1 and ECs 3 with interrupt 
operations and MC data transfer through an intcrcon- 60 signal transmissions on interrupt lines 8. 9 and 10. Mas- 
necting net of data and address buses. tcr computer input/output (MC I/O) unit U and scope 

Each of the multiple ECs performs a logic or arith- terminal 12 are conventional. The remaining element, 
metic function, being best characterized as individual manual control 13, will be specifically described herein- 
microprocessors. The MC, typifying a conventional after. 

scientific minicomputer, when actuated by the RTC 65 The schematic block by which the PDC structure 
while in the RUN mode, controls daU transfer among and operation are described herein are functionally 
all the functional elements, performs logic and arithme- labeled. Nevertheless, these blocks represent actual 
tic operations, and if called upon, performs the function apparatus cognizable by those practicing in the data 
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processing arts. The timing diagrams appearing in FIO. the EC function requiring the longest execution time. 
4 do not represent actual electrical wavefonns, but For instance, in one commercial microprocessor the • 
rather, are presented to show the time related interplay power function (Y^ requires an execution time more 
among the functional blocks in a qualitative sense. Ac- than 55 times greater than the simple multiplication of 
tuation of the interrupts is defined as occurring coinci* 5 Xx Y. 

dent with the waveform transition from a low to a high The other side of the balance allocating FT is the 
level. simulation layout, or stated otherwise, the distribution 

The MC as a whole, appearing in FIGS. 1 and 3 is of functions among the ECs, The PDC architecture 
broadly catagorized as a general purpose scientific digi- disclosed herein requires sequential processing and 
tal minicomputer having individual identifiable CPU, 10 transfer of data from EC to EC, with the MC CPU 
I/O and memory elements. It does differ in one minor designating the ECs to operation and the path of the 
respect Namely, the upper portion of the memory ad- data transfer. Since each EC must, in turn during its 
dresses, possibly 64^000 words, are reserved for conmiu- time segment, act upon the data processed by the previ- 
nication with the ECs. These registers are, nevertheless, ous EC during the previous EC time segment, and all 
still accessible to the MC CPU in the normal manner. In IS intermediate EC processed data must be available to the 
one version of the MC, it is further contemplated that integrator EC each FT, the maximum number of series 
the MC have the capacity to compile software as an connected non-integrator ECs between successive in 
adjunct to its PDC control functions. series integrator ECs is one less than the number of EC 

The ECs, 3, are primarily arithmetic processing units time segments in each FT. Again referring to FIG. 4. 
(APUs) with external interrupt capacity and various 20 the maximum number of non-integrator ECs between 
I/O registers. A simple APU would have a minimum of series integrator ECs is seven. In view of the above, it is 
two input registers, one output register and a command evident that the slower simulations, or ECs with faster 
register. Normally they are significantly more complex, execution time, provide more interconnection flexibility 
as shown by FIG. 2. Upon receiving the appropriate for the simulation programmer, 
interrupt signal from EC interrupt line 9 or 10. the EC 25 The RTC, when running, generates three chaimels of 
makes one pass of the selected program using the input output signals carried by lines 8, 9 and 10. The first line, 
register data, storing the results in the output registers. 8, goes to an interrupt on MC CPU 1. The second line, 
The program executed can either be one from a library 9, is directed to the ECs as an integrator interrupt. The 
of stored firmware programs, selected by the MC, or, in last line, 10, is the EC non-integrator interrupt As evi- 
a more exotic option, a program introduced by the MC 30 dent in FIG. 4, the MC CPU and non-integrator inter- 
for that particular pass. The pass internal to the EC is nipt time intervals are complementary in nature, i.e., 
synchronized to a separate, high frequency clock oper- there is no overlap of MC CPU and EC non-integrator 
ating at a rate commensurate with the command being operations. Each of these two time intervals is indepen- 
executed. As contemplated the EC I/O register data is dently adjustable, and the sum of the two substantially 
floating point format, though other data, such as pro- 35 equals the FT of the PDC system, 
gram selection commands, can be adequately handled As shown in FIG. 3, direct operator control over the 
with simple integers. RTC, with data entry access to the MC, is provided by 

The above-noted EC registers, in one form, are physi- manual control 13. Such external control is analogous to 
cally located with the individual ECs, for instance on •■ the operator panel of analog computers. As contem- 
the EC chip, but continue to have MC addresses and 40 plated, the simulation and RTC do not commence oper- 
access through the MC data bus concurrent with the ation until after a software compiler, which also runs on 
EC access. the MC, selects and stores the proper commands in each 

The RTC, as the name implies, is the pace setting EC command register, and generates the data transfer 
element for running the real-time simulation. It may or interconnect program to be run on the MC. 
generally be said that the RTC operates at a fixed rate 45 To adequately address the PDC architecture in a 
commensurate with the upper frequencies of the Simula- logical manner it is necessary to introduce optional 
tion, the execution times of the EC programs, and the design considerations, which to a substantial degree are 
largest number of non-integrating EC stages between cost effectiveness related. 

consecutive ECs functioning as integrators. Each as- The first of these design options is to load programs 
pect will be considered in turn at a later point 50 into the ECs from the MC, as opposed to having the 

The frame time (FT), as depicted on the plots in FIG. entire EC library stored on each EC for future selec- 
4, defines the time interval between integration updates. tion. This would be done during the PDC compilation 
Recognizing that the updates should preferably occur at and would replace the EC command mode selection 
a rate of no less than 10 times the highest frequency in introduced previously. A savings would be realized in 
the simulation, the nature of the simulation thereby fixes 55 terms of the memory required for each EC while also 
the maximum duration of the FT. As an example, if the eliminating the need for a command register in either 
simulation is that of an aircraft with responses up to 20 the EC or MC memory. On the negative side, compila- 
Hz, a FT of 5 msec, would permit integration updates at tion time would be increased somewhat. The second 
a rate 10 times greater, Le., 200 Hz. optional design feature offers tremendous potential for 

The minimum FT, on the other hand, is primarily a 60 the MC I/O; the option being an inclusion of the I/O 
tradeoff between the number and duration of the EC function in some or all of the ECs. Thereby, the EC 
time segments in each FT. The FT in FIG. 4 contains output could appear continuously to the user, though it 
eight such EC time segments with the remainder alio- would require that each EC library monitored have an 
cated to the MC and interrupt delays. As will become output program included. Outputs could be digital, with 
apparent from the example hereinafter, the major time 65 various formats, or D to A converted to permit direct 
consumers in each FT are the EC segments. attachments of actual hardware during simulation. A to 

On one side of the balance is minimum time that can D input data access into individual ECs is a corollary 
be allocated to each EC segment, this being bound by refinement Design option number three, the capability 
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to multiply each input by a constant, results in the sav- 
ings of many ECs when a complete simulation program 
is considered. The feature would aDow the PDC com- 
piler to load a constant at compilation time, which 
would multiply the result of the chosen program each 
EC time. If not needed, the PDC compiler would insert 
a value of one for the constant Experienced analog 
programmers can predate the savings in ECs as it 
would be analogous to saving an EC for each potenti- 
ometer in an analog simulation circuit The foiuth and 
last design option presented here is the inclusion of 
more than two input registers for some or all of the ECs, 
perhaps six inputs as shown in FIO. 2. This would result 
in a savings of ECs, particularly those performing sum- 
mation and multiplication functions. This does, never- 
theless, complicate the PDC compiler somewhat, so the 
tradeoff have to be considered in detail before imple- 
menting the technique into a PDC architecture. 
Sizing the PDC, i.e., determining the number of ECs 



only one stage transfer in the EC/MC register is desir- 
able. 

Under the PDC structure described above, integra- 
tion of a variable as a function of time is very straight- 
5 forward. Since synchronous operation is the modus 
operandi, the PDC is related directly to the integra- 
tion time interval and the integration algorithm be- 
comes minimal. For example, it is known in the art that 
three terms of the Taylor series expansion is a viable 
10 algorithm. 

la addition to the normal software for the MQ two 
special software items are needed, a PDC compiler and 
the individual program software for the EC library. The 
PDC compiler provides a means for translating a pro- 
gram written in a high level language, for instance a 
subset of Fortran with some special fimctions added, 
into the setup and initialization of the PDC The ma- 
chine language program generated controls the MC 
during the running of the simulation. Table 1 lists the 



15 



for such a computing system, is directiy dependent ^ principle functions of the PDC compiler as presently 



23 



30 



33 



upon the complexity of the simulation. For purposes of 
illustration, a representative count of functional de- 
ments in several six-degree-of-freedom analog simula- 
tion problems gave an average of 300 active compo- 
nents plus 300 potentiometers. If design option number 
three is employed, then 300 ECs is a good estimate for 
the minimum number and 1000 is a good maximum. 
Since it is assumed with this option that any of the EC 
hl)rary programs can be run on any EC, then the largest 
b1)rary program dictates the menoory size of the ECs. In 
further refmement, if design option number one is also 
employed, a good estimate is 4,096 words, with a maxi- 
mum of 8,192 words for EC memory. 

The amount of memory needed for the MC is dic- 
tated by both the complexity of the simulation and the 
extent of compilation necessary to implement the EC 
controls. Intuitively, 128K of 32 bit words appears de- 
sirable in the case of the six-degree-of-freedom simula- 
tion problems. In any case, the MC special memory ^ 
requirements are dir^y dependent upon the number 
of ECs used and the design option selected, particularly 
as to either one or four. 

Until now, '*words" as PDC computer elemaits have 
been mentioned, but word size has not been addressed 45 
as a design consideration. Based on analog accuracy, 
and fixed point arithmetic, 16 bit words for both the 
ECs and MC would appear to suffice. However, in the 
interest of eliminating scaling of variables floating point 
numbers are desired. This means that 24 bits for the so 
fraction and 8 bits for the sign and characteristic are a 
logical choice since the word size most commonly used 
b computers today is 32 bits, the next word size larger 
than 16 bits. 

As was suggested previously, the PDC should be 33 
fabricated as an intei^iited, hardware compatible de- 
vice. The logical approach to accomplish this is to first 
choose a master computer and then design and fabricate 
the ECs using the same technology. In the interest of 
speed, the MC data bus should be 32 bits wide. All 60 
arguments, parameters, data, etc., should be 32 bit float- 
ing point numbers to mgintflin resolution and eliminate 
scaling of variables. Of course, the MC and ECs should 
employ the same floating point format. The EC/MC 
registers located on the ECs, see FIO. 2, should be 
non-volatile with respect to the EC. Namely, input 
arguments to the EC should not be destroyed by output 
arguments from the EC. Also, in the interest of speed. 



contemplated. 



TABLE I 



EquBtkm Cncldcg 

Reduce oomplcx eqiuttoas into a number of equattoos 
or exprcsnoiu of simple form involving two variables 
and one aixthmetic operator. Tbe objective ta to 
divide tbe expressions so that each resultant exprea- 
uoa can be assigned to a single EC. 
Bo ok fc ptn g to allocate expressions to individual ECs. 
Selecting the commaDdi to designate the operating modes 
of each EC 

The commands are determined by Che arithmetic opera- 
tor in tbe simple form expression. 
Generating the MC RT program. 
The primary function of this program is to prescribe 
the transfer roate by which dau is routed from die 
ou^ut rcgistcn of each EC to the correct input 
regbten of other ECs. This b similar in function 
to the patchboard wiring of an analog computer. 
Oeoeraring a list of variables represented by the EC 



The individual software functions forming the EC 
library are foreseen to include any single or combina* 
tion of the functions listed in Table 2 below. 



TABLE 2 



Arithmetic Expressiotts 



subtraction 
mnltipUcation 
division 
exponentiation 

Logical Expressions 

not 

and 

or 

exclusive or 
Relational Expressions 
less than 
greater than 
equal to 
not eqnal to 



Functions 

All Fortran IV single Pre- 
dsioo non<complex func- 
tions 

Integration 

fixed interval 

Table Lookup 

one independent variable 



Actual operation of the PDC can be divided into 
three fairly distinct phases: program loading, program 
setup and initialization, and running the RT simulation 
program. During the loading phase, die high level lan- 
guage program is translated by the PDC compiler to 
63 allocate ECs, designate their operating commands, and 
generate a machine language program for controlling 
the MC during the running of the simulation program. 
Phase two. setup and initialization, consists of selecting 
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and setting starting values for the program variables, 
inserting data for table lookups, and selecting FT pa- 
rameters. It is desirable from an operator's viewpoint to 
be able to perform these functions via scope terminal 12. 

Each EC register 17, shown in FIG. 2, has a unique 3 
MC address accessible via buses 4 and 6 in addition to 
access from EC CPU 18. GeneraUy, all registers but 
those labeled OUTPUT can be written by the MC and 
read by the corresponding EC. The OUTPUT registers 
can be written by the EC and read by the MC >0 

Table 3 provides a fuiK:tional description of the activ- 
ities caused by actuating the various settings in manual 
control 13. It is worth noting that most arc analogous to 
contemporary analog computer operating modes in 
reasoning and sequence. 

TABLE 3 



CON. 

MODE FUNCTtONAL DESCRIPTION 



STOPS INTERRUPTS MC PROGRAM 
SETUP SO THAT IT WILL ACCEPT 

RTC DATA AND PARAMETER CHANGES 

STOPS INTERRUPTS MC PROGRAM SO 
RESET THAT IT LOADS DATA IN 

RTC PROPER MC/EC REGISTERS 23 

STARTS RTC RUNS MC TRANSFER 
RUN PROGRAM, EC COMMAND 

RTC PROGRAMS. AND INTEGRATOR EC 
PROGRAMS IN PROPER SEQUENCE 

STOPS INTERRUPTS MC PROGRAM SO 
HOLD THAT IT WILL ACCEPT ADDRESS 

RTC SELECTIONS TO MONFTOR EC 
OUTPUTS. 



For the ouster computer the priority, bigh to low, U: 
SETUP (Initiftltzation) 
RESET 
HOLD 
RUN. 



This means that SETUP mode can interrupt RESET, 
HOLD and RUN before their programs are completed; 
repeating in succession as to the remaining modes. Note 40 
also, if a MC program is interrupted, the interrupt sys- 
tem does not cause execution of the remaining program 
at a later time. These interrupt signals are routed inter- 
nally to the EC in accordance with the specific firm- 
ware program designated by the data in the COM- 45 
MAND register. 

Actual operation of the PDC commences after the 
SETUP mode has been fully executed by entry and 
compilation of the real-time simulation program. Upon 
actuation of the RUN mode, EC programs in the read so 
only memory (ROM) of each EC, as selected by their 
COMMAND registers, are executed upon the appropri- 
ate interrupt signal from the RTC. During the MC 
segment of the RTC, data is transferred from the EC 
output registers to the appropriate EC INPUT regis- 55 
ters, being substantially analogous to the interconnects 
on analog computer patchboards. Though this opera- 
tion is performed sequentially, it is one of the faster of 
digital computer operations. All ECs, except the inte- 
grators, execute their programs simultaneously or in 60 
parallel upon receipt of the interrupt signal, repeating as 
depicted in FIG. 4 over the duration of the FT. At the 
end of each FT all integrator ECs are updated in paral- 
lel, following which the RTC commences a new FT. 
The process repeats cyclically as long as a higher prior- 63 
ity interrupt is not initiated from the manual control 13. 
Changes in mode conform to the selection outlined in 
Table 3. 



.025 
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To better illustrate the method by which a PDC 
structure is defined, consider the following illustrative 
example. In this case the simulation to be performed is 
that of a six-degree-of-freedoro aircraft or missile. Table 
4 outlines the elements needed. 

TABLE 4 

Namber of 



Type of EC 


Kumber 
of ECs 


MC/EC 
Registcis 


Required 
Addresses 


Gexteni] purpose (GP) 


170 


7 


1.190 


GP with A/D 




7 


4Z0 


GP with D/A 


60 


7 


420 


Iotegr8tx>rs 


190 


7 


1,030 


Function Generator! 


60 


1004 


60,240 






TOTAL 


63320 



As suggested previously, the MC memory should con- 
tain 128 K of addressable 32 bit words. From Table 4 it 
is evident that 64 K of the memory need be reserved for 
MC/EC registers while the remaining 64 K is available 
as conventional memory. 

In fairly conventional manner the timing signals 
would preferably begin with a master clock of high 
frequency, at the cycle rate of the MC and ECs, cycle 
counted to generate the RTC signals. 

The geneal purpose ECs in Table 4 have register 
configuration 21 of FIG. 2 and perform the functional 
computation: 

OUTPUT=aA.B, CD). 

where function f is the one from Table 5 which corre- 
sponds to the COMMAND register data. A, B, C. D 
and E are names given to the INPUTs. They can repre- 
sent problem variables and change during RUN mode 
or represent constant coefficients, very similar to poten- 
tiometers in an analog circuit, that are prescribed during 
RESET mode and remain unchanged during RUN. If 
no multiplying is needed, or if only one input is desired, 



the coefTicient is set to one or zero. 

TABLE 5 


COMMAND 


OUTPUT FUNCTION 
(FLOATINO POINT) 


ADD 


A X B + C X D 


SUB 


A X B - C X D 


MUL 


A X B X C X D 


DIV 


A X B 4- C X D 


SORT 


lA X B + C X D]l 


SIN 


SIN [A X B + C X D] 


COS 


COS [A X B + C X DJ 


TAN 


TAN [A X B + C X DJ 


ASIN 


ARCSIN (A X B + C X DJ 


ATAN 


ARCTAN [A X B + C X Dl 


LOG 


IjOGio(A X B + C X di 


LN 


LOG, [A X B + C X Dl 


EXP 


JA X B + C X Dl 


PWR 


(A X B)tCxXI> 


ASW 


larger [A X B] or IC X Dl 


LSW 


A X B E > 0 




C X D E < 0 



If the EC is to accept analog data for one of the in- 
puts, for instance D, it enters on line 22 in the analog 
state and is converted internally to the digital equiva- 
lent. The function expression for the ADD command is 
written as: 



OUTPUT = A X B + C X AD. 
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Conversion of the output to analog is performed in a 
similar manner with the analog signal appearing on line 
23. 

The ECs functioning as integrators perform a calcu- 
lation which approximates the analytic expression: 



/ XdL 



The algorithm used is: 

Ynrw°fTX YcU^FTniX^-Xotd) 

Xin is the data entering the EC integrator stage, while 
COEFF is a scaling multiplier analogous to an input 
potentiometer in an analog system. Together they pro- 



MC address buses, the latter being for the purpose of 
storing the data directly from the MC. 

Implementing the PDC requires the consideration of 
three crucial operating constraints. 
5 (1) Can the EC programs be completed in the time 
segments allocated by the RTC? 

(2) Is there is a computational lag or error from one 
EC to another due to the hardware configuration? 

(3) Is the overall computer speed sufficient for real- 
10 time simulation? 

Although it is not shown explicitly in the figures 
depicting the block diagram of an EC. all APUs of such 
ECs have feedback signals indicating that the execution 
of the previously entered command is completed. These 
IS signals are used in a hardware check system to initiate 
an alarm if a RTC interrupt initiates a new command 
prior to the completion of previous commands by all 
ECs. 

The possibility for computational lag. the second 



vide added flexibmty to each EC integrator. The pa- ^ ^^^^^ joes exist based on the PDC architecture 



rameters FT and FT/2, as well as the initial values of X 
and Y, enter their respective registers when the RESET 
mode is actuated. Initially, the registers are set such 
that: 



Ym 



' yinitiai 



disclosed here. However, the provision for avoiding the 
lag is also inherent in the architecture and it would be 
the responsibility of the programmer to insure that it is 
avoided. A sophisticated PDC compiler can also be 
2^ used to test for this possibility. The rule for avoiding 
computational lag due to data transfer between ECs. is 
to count the largest string of other ECs between any 
two successive integrator ECs within the program, and 
choose a cycle time which exceeds that count by at least 
^ one. It is anticipated that this count will seldom exceed 
ten. and based on the timing information which follows, 
will permit operation without error. 

Timing estimates for the particular example under 
consideration are shown in Table 6. and indicate that no 
This transfer step means that after and \new are 35 problem exists for the real-time simulation as it is de- 
— r^ii^.,^-- .v-. »i — - „i — A»*» Y practicing in the art recognize 



Xotd-Xafw^'XinHfai, 

whereupon they are updated following the sequence: 



40 



calculated following the above algorithm, data Xnfwand 
Ynm^are stored in registers Xo&/and Yo&/for the succeed- 
ing FT calculation. The EC register configuration for 
integration is shown by reference numeral 24 in FIG. 2. 

Single variable functions can also be readily accom- 
modated in the PDC architecture. Like most computer 
circuit designs, speed and cost are prime factors that 
influence the end product, and this is no exception. It is 
well known that, if sufficient data is stored, no interpo- 
lation between data points is necessary. Obviously, this 
is a tradeoff between memory and computation time. 
Since the cost of random access memory (RAM) is low. 
the large storage alternative to single variable function 
generation b appealing. Based on experience. 1000 data 
points used to represent aerodynamic functions em- 
ployed in aircraft or missile simulations would appear to 
allow storage of data sufficiently close to omit interpo- 
lation between data points. With such an approach. 
Xmaxf XmiA and the intermediate X value data are intro- 
duced into storage during the RESET mode. During 
the RUN mode values of the function, f(X). are deter- 
mined by the data stored in registers. The registers are 
addressed by integers such as IX, here IX is the integer 
value for the floating point expression: 



that improvements in speeds of microprocessor devices 
arc literally forthcoming daily, so the estimates in Table 
6 are comparatively conservative. 

TABLE 6 

Microprocessor cycle clock frequency is lOMHz ' 
1 ckx:k cycle » 10~^ seconds 

Typical EC instruction (add. sub. mul, div) execution cycles 
200 clock cycles or 

200 X 10-' « 2 X lO-^scc/instnictioo 
Worst cose EC mstructkm (POWR) execution cycles 
9SO0 clock cycles or 
10* X 10"' sec/Instruction 

Interrupt service time for each interrupt (assuming the wave- 
form in FIO. 4 applies) 
50 8 clock cycles or 

8 X I0-' sec/interrupt 

Max desirable FT for real-time simulation having frequencies 
up to ^ Hz. 
5 X IQ-^sec/FT 
Number of typical instructtons/FT 



55 



5^10-3sec/FT ^ ^50 instr/FT 



2 X 10-' sec/instr 



Number of worst case instrvctioos/FT 



60 



5^10-3sec/Fr ^3i„^,^ 



fX B {Xij, — Xmia) when 



I ~ Xifi 



6S 



10-^ sec/instr 

NOTE: 

This is KkkJoi used A mighi even be drepped from inuracf ion lei. 
Max time used by interrvpt/FT (assuming FIG. 4 applies) 
250int/FT X 8 X 10" ' sec/btcnupt = 2.0 X lO-^sec/FTc 



As shown in reference numeral 26 PIG. 2, the data in 
the registers can be addressed by either the EC or the 



2.0 X 10-^ 
5 X 10-^ 



X 100 = 4.0% of FT 
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Integr&tor update time 

10^ cyclq/mtcgration X tO~^ >cc/cyclc « 10" 



From the calculations presented in Table 6, it is clear 
that aU the constraints can be satisfied without undue 
difTtculty if moderate care is exercised in the distribu- 
tion of non-integrator ECs between successive integra- 
tor ECs. This is the providence of the programmer, but 
may, as previously noted, be checked by the MC during 
compilation. 

To further expand upon the illustrative example of 
this embodiment, though with lesser complexity, con* 
sider the simulation undertaken to be that of a point 
mass missile attempting to hit a point mass target, where 
the target is traveling at a constant velocity and turn 
rate. For purposes of this limited- example, motion is 
constrained to the horizontal plane with proportional 
navigation as the guidance law. This example demon- 
strates the FDC programs, the sequence of events nec- 
essary to solve the problem, and a direct comparison of 
EC and analog computer elements. The symbols, the 
parameters they represent, and the corresponding units 
used in the example are defmed in Table 7. 

TABLE 7 



Computer 
Equatioa Program 



Parameter 



Symbols 


Symbols 


Represented 


Units 


R 


R 


range, missile to target 


feet 


X 


LAMB 


line to sight angle to X 
reference 


degrees 


Jm 


GAMM 


angle of missile velocity to X 
refcretice 


degrees 


It 


OAMT 


angle of target velocity to X 
reference 


degrees 


v« 


VM 


missile velocity 


ft/sec 


V/ 


VT 


target velocity 


n/sec 




XM 


position of missile on X axis 


feet 


v„ 


YM 


position of missile on Y axis 


feet 


X, 


XT 


position of target on X axb 


feet 


Y, 


YT 


position of target on Y axis 


feet 


t 




time 


seconds 


K 


K 


gain constant for steering 




o 


O 


subscript representing initial 
value 

suffix representiog initial value 






FT 


frame time 


seconds 




RTI 


RTC interval to MC 






RT2 


RTC interval to oon'tntegratof 
ECs 





Table 8 contains the analytical equations defining the 
example simulation. Since the equations are well known 
to those performing simulations in the aerospace field, 
and are extensively described in publications, their ori- 
gins will not be described. 

TABLE 8 

R « V|Oos(\ -yd- V;„cos{X - y„) 

R = / R dt + R<, 
o 

X = [-V/on(X - y,) + V„,sin(X - 'y„)ll/R 

X- /*Xdi + X<, 
o 
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IS 



20 



ym « / Tmdt + 
o 

X«« V^cosX* 
Y„ - V„ sin X„ 



Y„ = / Y«dt + Y„ 



X » V,co» yt 
Y, = V,siny, 

X, = / X/li + Xtt 



Y, « / Y4t + Yft, 
o 

FIG. 5 is a schematic diagram of an analog simulation 
25 corresponding to the set of equations in Table 8. As 
shown, the analog diagram has been simpliHed to the 
extent that no scaling is done and sin/cos generators are 
asstmied capable of processing either positive or nega- 
tive inputs. Also appearing in FIG. 5 are overlaid 
30 blocks, defined by d^hed lines, enclosing the various 
analog system functional elements. Each block repre- 
sents an EC performing the functions enclosed. Again, 
since the simulation is comparatively simple and well 
known by those practicing in the art, FIG. 5 merits no 
35 further elaboration. 

FIG. 6 is the same simulation diagrammatically pres- 
ented in PDC block diagram structure. The blocks 
represent individual ECs while the dashed intercon- 
necting lines define the MC transfer program, i.e., the 
40 flow of data between ECs under the control of the MC. 
For purposes of comparison, the blocks in FIGS. 5 and 
6 are correspondingly numbered in parenthesis. To 
relate the ECs to the computer program described here- 
inafter, the MC addresses for the registers follow a 
43 format of EC nni. The nn represents a two digit sequen- 
tial EC number and i represents a MC/EC register 
associated with EC nn. The registers conform to those 
appearing in FIG. 2, reference numerals 21, 24 and 26, 
with those having seven registers numbered 1 through 7 
50 proceeding downward. 

A full PDC program for the example simulation is 
provided in Table 9. The numerical values selected are 
considered to be nominal. Since the equations, symbols, 
and interconnected block diagrams have been intro- 
S5 duced, the program itself is substantially self-explana- 
tory with minor exceptions. The program under 
SETUP IS inserted when the PDC is in that mode. The 
RESET program is entered during its corresponding 
mode, though entry would likely be directly from inter- 
60 nal MC storage programmed prior to initiating RESET. 

TABLE 9 

RUN 

SETUP PROGRAM RESET PROGRAM PROGRAM 



63 



OAMT = 3.1416/2. EC0I6 = ADD ECOII = ECI07 

RO - 10000. EC026 = ADD EC02I = ECIQ? 

LAMBO»a EC036 = COS £0023 » ECU? 

OAMMO « a EC046 « ADD EC031 = EC027 
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SETUP PROGRAM RESET PROGRAM 



RUN 
PROGRAM 



XMO « a 


EC066 M SIN 


£0041 


M EC1B7 


YMO » a 


EO076 a SIN 


EO043 


a £0037 


XTO » loooa 


EO086 — ADD 


EOOSl 


« EO047 


YTO m a 


EO096 M DtV 


EO061 


BO017 


FT « .001 


ECl6fi a COS 


E0071 


« EOQ27 


RT1 « .0001 


EC176 * SIN 


EOOSl 


M EOD67 


RT2 » .0001 


EC1S6 « COS 


EO083 


M EO077 


VM - 9oa 


EO012 « 1^ 


EO091 


EO087 


VT = 700. 


EO013 « -OAMT 


EO093 


n'ECQ37 


K - 4 


EO014 • 1.0 


EClOl 


« £0097 


BSD 


EOQ22 1.0 


EClll 


» EO097 




EO024 ^ -t.0 


EC12I 


« EC167 




EC032 - 1.0 


EC131 


■» EC177 




EO033 B a 


EC161 


w EC117 




£0042 « +VT 


£0171 


» E0117 




ECOa - -VM 


£0181 


- £0017 




EO033 - FT 


END 






EOOM » FT/2. 








EO03S » 1.0 








EO0S6 - RO 








EO037 RO 








EO062 » IjO 








EO063 - a 








Eoon » 1.0 








£0073 - a 








EO082 » -VT 








EO084 « VM 








EO092 » 1.0 








EO094 - 1.0 








EC103 » FT 








£0014 » FT/2. 








ECtOS - 1.0 








EC 106 s LAMBO 








ECI07 t. LAMBO 








EC113 » FT 








EC114 » FT/2. 








ECtl3 = 1.0 








BCU6 a OAMMO 








ECU7 «- OAMMO 








EC123 » FT 








EC124 = PT/1 








EC123 *a 1.0 








EC 1 26 H XMO 








EC 127 - XMO 








EC133 — FT 








EC 134 ■> FT/1 








EC135 « li) 








ECI36 « YMO 








EC137 a YMO 








£0141 « 
















£0143 B FT 








EC144 « FT/2. 








BC14S K 1.0 








EC 146 M XTO 








ECI47 - XTO 








EClSl » 








VT«SIN(OAMT) 




i 




ECI53 - FT 








ECI54 - FT/2. 








EC135 » 1.0 








E0136 « YTO 








EC157 ^ YTO 








EC162 - 1.0 




i 




EC163 a a 








Ecin « li) 








EC173 - 0. 








EC1S2 - 1.0 








EC183 «* a 








END 




* 



One skilled in the art clearly recognizes the plethora 
of attributes distinguishing the PDC architecture dis- 
closed herein from the conventional analog computer. 
As to physical characteristicSp the conventional patch- 
cord interconnect system and the high maintenance 
costs of the analog version have been Mimmfif^ Qood 
accuracy and full repeatability are attained from the 
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digital format and its resolution, while simultaneously 
eliminating the scaling of variables burden present in 
analog systems. As with all digital systems, program- 
ming flexibility and large data storage capability are 
inherent characteristics. Direct interface with analog 
hardware is not sacrificed since A/D and D/A conver- 
aon devices are readily available. And finaUy, present 
estimates place the purchase cost in the range of half 
that of an analog system while having the capability to 
support over 90% of the typical simulations undiertaken 
on the comparable system. 
I claim: 

1. A digital data processing apparatus for executing 
real-time analog format functions comprising: 

(a) a plurality of elemental computers, each of which is 
a complete digital data processor having its own 
central processing unit, memory, input/output 
means, cycle clock, external interrupt capacity, and 
stored program; 

(b) a master controller with a master cycle clock and an 
external interrupt capacity for regulating the transfer 
of data between elemental computers in conformity 
to an ordered sequence; 

(c) an MC memory; 

(d) master bus means comprising a data bus and an 
address bus, both coupled to the master controller, to 
the MC memory, and to the input/output means of 
the elemental computers; 

(e) a real-time clock operating at a rate commensurate 
with the upper frequencies of a simulation but at a 
substantially lower frequency than said cycle clocks, 
an M C interrupt line coupling said real-time clock to 
the master controller, and two MC interrupt lines 
which are an integrator interrupt line and a non-inte- 
grator interrupt line, with each elemental computer 
coupled to said real-time clock via at least one of said 
two EC interrupt lines; wherein said real-time clock 
generates timing comprising repetitive frames with 
alternating EC time slots and MC time slots in each 
frame, the EC time slots being divided into a plurality 
of non-integration time slots and one integration time 
slot each frame, the number of EC time slots per 
frame being determined by the ""^if n^Tm number of 
interchanges of data required by the program for 
non-integrating functions between integratmg func- 
tions; 

(0 wherein for any given simulation program some of 
the elemental computers are assigned non-integrating 
functions and are enabled via said non-integrator 
interrupt line during non-integration time slots, and 
other of the elemental computers are assigned inte- 
grating functions and are enabled via said integration 
interrupt line during integration time slots, each said 
function corresponding to the function of one com- 
puting element of an analog computer, each elemen- 
tal computer making one pass of its program to per- 
form its assigned function digitally upon receiving an 
interrupt agnal; 

(g) the master controller being enabled via said MC 
interrupt line diiring MC time slots to transfer data 
via said master bus means among the elemental com* 
puters, so as to be means functionally equivalent to 
patchcord interconnections of computing elements of 
an analog computer, and 

(h) wherein said input/output means comprises registers 
addressable from either the master controller via said 
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master bus means or from the respective elemental lion of the remaining program at a later tunc; and the 

computers, setup mode has priority to interrupt the reset mode. 

2. The apparatus as recited in claim 1 wherein said 5. The apparatus as recited in claim 4, further includ- 
mastcr controller is a general purpose scientific digital ing a hold mode which is means operative to interrupt 
minicomputer having its individual identifiable master 5 the run mode so that the program will accept address 
central processor, master input/output means, ialemal selections to monitor the elemental computer outputs, 
interrupt capacity, and memory elements included in the hold mode being interruptable by either the reset or 
said MC memory, differiiig in one respect in that a the setup mode. 

portion of the memory addresses are reserved for com- 6. The apparatus as recited in claims 1, 2, 3, 4 or 5 

munication with said elemental computers, wherein at 10 wherein said input/output means of some of the elemen- 

least some of the elemental computers contain internal tal computers includes four input registers A, B, C, D. 

thereto a library of processing functions selectable by wherein registers B and D are used to store constants 

the master controller via a command register forming which remain fixed in value as long as the real-time 

part of said input/output means, and wherein said mas- clock is running, with the elemental computer ftinction- 

ter bus means is the only connection of each elemental 13 ing to multiply the contents of registers A and B and 

computer to the master controller and to other elemen- q d to form two inputs for the principal 

tal computers for data and instruction transfer, the function that elemental computer is programmed to 

transfer being only to. from or via the master controller. perform, to thereby greatly reduce the number of cle- 

3. The apparatus as recited in claun 2, wherein the mental computers required. 

master controller is able to run general purpose digital 20 -j^^ apparatus as recited in claim 6, where said 

programs and also has the capacity to compile software input/output means for integration elemental comput- 

as an adjunct to its control fimctions. ^^s includes registers for old and new values of depen- 

4. The apparatus as recited in claim 2, wherein there ^^^^ independent variables, one register for a coeffi- 
is means for providing that operation is divided mto cjent and at least one register for a frame timing paramc- 
three fairly distinct phases designated respectively as a 25 

setup mode, a reset mode, and a run oKxle; a manual j apparatus as recited in claim 7. wherein the 

control analogous to the operator panel of an analog input/output means of some of the elemental computers, 

computer for selection of the mode and contiolhng jndude a large number of registers for storing values of 

insertion of data; fiinction of one variable. 

the setup mode bemg means operaUve for program 30 ^ apparatus as recited in claim 8, wherein the 

loading to allocate the elemental computers, desig- t/^^j j ^^ans of at least some of the elemental 

nate their operatmg commands, and generate the ^^ ^^^ -^^^^^ eternal mput/output means inde- 

program for controUmg the master pendent of said bus means, so that the elemental com- 

the run mode, said real-ume clock bemg stopped ^t^^ output may appear continuously to a user, inlud- 

dunng this mode; ^ ing digital-to-analog conversion for some of the out- 

thc reset mode bemg means °P^^^«;°^ ?f^"P Jts, L analog-tcAligital conversion for some of the 

initialization, comprising sclectmg and settmg start- puia, oiui » « 

ing value* for the program variables, inserting data "^^^^'^ ^ 1, 2 or 4 

2d" ^jsrci^^^^^^^^ « 

7^ . . elemental computers mcludes external data transfer 

thrr^oSS m<^;.erative with said real-time m«u. Independent of said ^ ^^J^^.. . 

clock running fo run thrnon-integrator elemental The apparatos <|s recited m chum 10, farther m- 

^mpuws simultaneously during EC time slots, the eluding analog-t«><i.g.tal conveiswn mean, assocated 
^e^^ntroller daU tLsfer Jrogram during MC « with saW extemd data transfer mea^ ^^^^ 

time slots, and during the integration time slot to mput of data to some of the elemental «>«P« «?;^ 

update in parallel all integrator elementary comput- i2. The apparatus as recited m clami 11. further m- 

g{~ eluding digital-to-analog converter means associated 

with interrupt means capabUity wherein either the setup with said external data transfer means to permit direct 
or the reset mode has priority to interrupt the run SO attachment of actual^hardware^dunng smiulaUon. 

mode, and the interrupt system does not cause execu- * 
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